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SUMMARY 

A probabilistic approach to composite micromechanics is developed that 
simulates uncertainties in unidirectional fiber composite properties. These 
methods are in the form of computational procedures using Monte Carlo simula- 
tion. The variables for which uncertainties are accounted include constituent 
and void volume ratios, constituent elastic properties and strengths, and 
fiber misalignment (primitive variables). A graphite/epoxy unidirectional com 
posite (ply) is studied to incorporate fiber composite material properties 
uncertainties at the micro level. Probabilistic composite micromechanics pro- 
vides extensive information which formally relates ply phenomenological behav- 
ior to a large number of complex and interacting uncertainties at the 
constituents 1 evel . 


INTRODUCTION 

Properties of composite laminates depend on the properties of the constit 
uent materials, their distribution, and orientation. Laminates are composed 
of layers of unidirectional ly reinforced plies (lamina). The ply is typically 
considered the basic unit of material in composite mechanics. The branch of 
composite mechanics that predicts ply material properties based on the proper- 
ties, volume, and orientation of its constituents is known as composite 
micromechanics, and frequently incorporates the traditional Mechanics of 
Materials assumptions. The desired laminate is created by stacking piles in 
specific directions. The integration of ply properties to yield laminate prop 
erties is called laminate theory. Laminate variables such as ply orientation 
and stacking sequence can be tailored to yield a laminate with the desired 
material properties. Thus, the laminated composite is a suitable material for 
component design. 

Analysis of fiber composite structures is currently performed using a 
variety of computer codes. Complete mechanical, thermal, and hygral proper- 
ties are calculated, and can be used to compute respective structural 
responses. Advanced failure criteria are used to calculate composite 
strengths. Environmental effects are also quantified. The usefulness of 
these codes has been demonstrated by comparison with experimental data and 
finite element results, (ref. 1) 
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The analytical capability of these codes is limited by the deterministic 
nature of the computations, while corresponding measured data exhibit consider- 
able scatter. Specifically, fixed values for constituent material properties, 
fabrication process variables (i.e., constituent volume ratios), and internal 
geometry (primitive variables) are used as input. However, random variations 
in these parameters are not only expected, but easily observed experimentally. 

The analysis of composite structures requires reliable predictive models 
for material properties and strengths. However, the prediction efforts have 
been complicated by inherent scatter in experimental data. Since uncertainties 
in the constituent properties, fabrication variables, and internal geometry 
would lead to uncertainties in the measured composite properties, the question 
arises: 

How much of the "statistical" scatter of experimentally observed composite 
properties can be explained by reasonable statistical distribution of 
input parameters in composite micromechanics and laminate theory predic- 
tive models? 

In order to answer this question, a study was conducted to develop a pro- 
babilistic approach to composite micromechanics. The objective of the present 
paper is to describe this approach and present typical results. The computa- 
tional simulation is preformed using ply substructuring with an existing com- 
puter code (ref. 2) for composite mechanics and in conjunction with Monte Carlo 
simulation. The randomness in the primitive variables is selected from 
anticipated respective probabilistic distributions. These distributions have 
the same mean as the deterministic case. In this context, the deterministic 
micromechanics will predict only the mean while probabilistic will predict the 
mean and respective scatter. 


SYMBOLS 

E normal modulus - subscripts define type and direction 
FVR fiber volume ratio 

G shear modulus - subscripts define type and direction 

S fracture stress (strength) - subscripts define type and direction 

VVR void volume ratio 
Greek 

a thermal expansion coefficient - subscripts define type and direction; 
also shaped coefficient the Weibull distributions 

(5 mean value in the Weibull distribution 

y special probabilistic distribution 

k shape coefficient in the r distribution 
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X 


mean value in the r distribution 


M mean value in the normal distribution 

v Poisson's ratio - subscripts define type and direction 

<7 standard deviation 

Subscripts 

c compression 

f fiber property 

51 ply property 

M matrix property 

S shear 

T tension 

1 fiber direction and surface normal to this direction 

2 transverse fiber direction and surface normal to this direction 

3 thickness direction and surface normal to this direction 


COMPUTATIONAL SIMULATION FUNDAMENTALS 
Analysis Model - Basic Unit 

The models commonly used in characterizing fiber composites structural 
behavior are based on the calculation of properties of the basic unit of an 
orthotropic ply. The layup geometry is then used in laminate equations to cal- 
culate composite properties (fig. 1). In the present work, however, the basic 
unit is taken as a sub-ply (ply substructuring), which consists of only one 
fiber-matrix level in the material (fig. 2(a)). Micromechanics theory is used 
to calculate the properties of the assumed orthotropic sub-ply, each with ran- 
domly distributed fabrication variables (fiber volume ratio and fiber misalign- 
ment) and material properties. Uncertainties in fiber directions, due to 
possible misalignment within the ply, are then used in the laminate equations 
to calculate ply properties. The substructuring of the ply represents a novel 
attempt at characterization of fiber composite material properties based on 
probabilistically distributed constituent properties, individual fiber misal- 
ignment (fig. 2(b)) and fabrication process (primitive) variables. 

This formulation is particularly well suited to probabilistic description 
of fiber composite material properties. The micromechanics and laminate theory 
equations can be used to calculate ply properties at any number of points in a 
ply. The approach provides a computationally rational procedure for evaluat- 
ing the scatter in composite material properties because it evaluates ply 
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behavior as the result of a series of random occurences (uncertainties In the 
primitive variables) which occur at the intraply or micromechanics level. 


Composite Mechanics 

The simulation Is performed by considering the ply as an assembly (equiva- 
lent laminate) of 15 subplies. The composite mechanics used in the simulation 
Is that embodied in the Integrated Composite Analyzer (ICAN) (ref. 2), which 
Is a computer program for comprehensive linear analysis of multilayered fiber 
composite structures. The program contains the essential features required to 
effectively design structural components made from fiber composites. It now 
represents the culmination of research conducted since the early 1970's, at 
the National Aeronautics and Space Administration (NASA) Lewis Research Cen- 
ter, to develop and code reliable composite mechanics theories. Through this 
simulation the following are Included and can be evaluated probabilistically: 

(1) Conventional laminate analysis 

(2) Intraply and Interply hybrid composites 

(3) Hygral, thermal, mechanical properties and response 

(4) Ply stress-strain influence coefficients 

(5) Microstresses and microstress influence coefficients 

(6) Stress concentration factors around a circular hole 

(7) Predictions of delamination locations around a circular hole 

(8) Poisson's ratio mismatch details near a straight free edge 

(9) Free edge Interlaminar stresses 

(10) Laminate failure stresses 

(11) Normal and transverse shear stresses 

(12) Explicit specification of matrix-rich Interply layers 

(13) Finite element material cards for general purpose finite element 
structural analysis 


Probabilistic Simulation - Monte Carlo Methods 

Complicated stochastic (probabilistic) processes can be simulated by a 
variety of numerical methods generally referred to as Monte Carlo methods 
(ref. 3). The term refers to that branch of mathematics concerned with numeri- 
cal experiments on random numbers. Since the advent of high speed computers, 
they have found extensive use In most fields of science and engineering, In 
analyzing many physical processes of a probabilistic nature, or where physical 
experimentation is not feasible. In general, they can be economically used to 
achieve a level of confidence of greater than 90 percent. 
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A Monte Carlo simulation refers to the procedure of randomly assigning a 
value to an independent random variable in a chosen model, and observing the 
dependent variable at the conclusion of the process being modeled. A Monte 
Carlo simulation is composed of n such independent experiments. When n is 
sufficiently large, the observations will yield a statistically meaningful 
estimate of the model dependent variable. 

The form of Monte Carlo used in the present investigation is as follows: 

(1) Define the system model by assuming: 

(a) Probability distributions for the scatter in all independent 
(primitive) variables. 

(b) The equations used to model the composite thermal and mechanical 
behavior properly describe the physics. 

(2) Use the random sampling techniques to select input values of the inde- 
pendent variables from their assumed probabilistic distributions 
(1(a) above). 

(3) Calculate dependent (output) variables using the equation. (Kb) 
above) . 

(4) Estimate regression parameters (goodness of fit tests) for the assumed 
model . 

(5) Replicate the procedure, each time with a new set of input values, 
using the procedure in (2) above. 

(6) Use appropriate statistical methods to calculate properties of the 
dependent variable distributions such as means and standard varia- 
tions (results of this step are not included in this paper). 


Computational Simulation Procedure 

'o perform the computational simulation, a computer code was developed to 
ICAN and an available statistical analysis code (ref. 4). The logic 

m for this code is shown in figure 3. The steps are as follows: 

I) Select values for the primitive variables for each subply from their 

respective assumed probabilistic distributions: 

(a) Normal with selected variations - constituent elastic properties 
and fiber volume ratio, and fiber misalignment 

(b) Weibull - constituent strengths 

(c) Gamma - for void volume ratio 

Fifteen different sets are generated (one for each subply) where the 
means and variation ranges selected for the fourteen primitive random varia- 
bles are those typical for AS-graphite fiber/epoxy composites as listed in 
table 1. The deterministic case is also shown. Note that the deterministic 
case values are the same as the means for the other two cases. 
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(2) Enter these values as inputs into ICAN. 

(3) Run ICAN and retrieve and store ICAN output for desired ply 
properties . 

(4) Repeat process n-times where n is sufficiently large to provide 
data with an acceptable level of confidence. For the 

results presented here n equal 50 was considered adequate. 

(5) Graphically stored output for probability density and cumulative 
probability distributions. 

Additional details and rationale are described in reference 5. 


RESULTS AND DISCUSSION 

Results obtained from probabilistic composite micromechanics are volumi- 
nous. These results can be presented in several alternative ways. A meaning- 
ful way to present them is in terms of (1) histograms (frequency or probability 
of density) and (2) cumulative distribution. The histograms indicate the 
range of the scatter while the cumulative distribution indicates the probabil- 
ity of occurence. The graphical results presented for the ply properties are 
(1) moduli (longitudinal, transverse, shear and Poisson's ratios, (2) thermal 
expansion coefficients (longitudinal and transverse), and (3) strengths (longi- 
tudinal tension and compression, transverse tension and compression and intral- 
aminar shear). For comparison purposes the deterministic values for these 
same properties are listed in table II where the inputs for the deterministic 
values are from Case 1, table I. 


Modul i 

Longitudinal modulus (E on) - The frequency and cumulative distribution 
diagrams for the longitudinal modulus are shown in figure 4. Case 2 and Case 
3 refer to probabilistic inputs in table I. As would be expected, the greater 
the scatter in the primitive variables, the greater the scatter in the ply 
modulus Case 2 versus Case 3. For example, the scatter in the longitudinal 
modulus for Case 2 is from about 13.5 to 18.5 mpsi while that for Case 3 is 
from about 12.5 to 19.5 mpsi. The estimated means for the respective cases, 
taken at 50 percent probability of occurence, are about 15.5 and 14.8 mpsi. 

Both of these are less than the deterministic case of 15.8 mpsi (table II). 

A typical experimental value of modulus for this composite system is 
about 16 mpsi which is well within the range of the scatter for both cases. 
Assuming the scatter for the primitive variables is reasonable, the results 
show that experimental values for the longitudinal modulus as low as about 
12.5 mpsi and as high as 19.5 mpsi are probable under present practice. 

Transverse modulus (E p 77) - The frequency and cumulative distribution dia- 
grams for E&22 are shown in figure 5. Focusing on the diagrams for Case 3, 
it is seen that the scatter for E^z is from the 1 to 1.3 mpsi with a mean 
value of 1.1 mpsi. This mean value is about the same as the 1.06 mpsi for the 
deterministic case. A typical experimental value is about 1.5 mpsi, which is 
greater than the highest value of 1.3 mpsi in the scatter range. One possible 
explanation is that the mean value used as input for the fiber 
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transverse modulus is lower than its actual value. This is readily corrected 
by a direct shift of the mean in the input value for Eg,22 » or a mean ra tio 
shift of (1. 5/1.1) times the imput value for Eq,22- 

Assuming the scatter used for the primitive variables is reasonable, 
then, the anticipated scatter in measured values for Eg,22 is about 20 per- 
cent. Variations greater than 20 percent are likely to arise from factors 
which were not considered in this investigation and which are not that obvi- 
ous. This illustrates another important benefit of probabilistic composite 
mi cromechani cs . 

Shear modulus (G p.i?) - The frequency and cumulative distribution diagrams 
for G a i2 are shown in figure 6. One significant observation from these dia- 
grams is that the means for the two cases differ by about 0.2 mpsi or about 
30 percent. The scatter assumed for the primitive variables significantly 
affects the anticipated mean value of the ply shear modulus. 

Looking at the diagrams for Case 3, we see that the scatter for Gq,i 2 
varies from about 0.6 to 1.4 mpsi with a mean of about 0.81 mpsi. The deter- 
ministic value is 0.52 mpsi (table II). The mean value is about the same as 
the lowest value in the scatter range. Three observations follow: 

(1) additional sampling points need to be included in plotting the diagrams, 

(2) the assumed scatter in the primitive variables (especially fiber misalign- 
ment) tends to increase the ply shear modulus relative to the deterministic 
case, and (3) a large scatter is probable in measured values for Gq, 12. This 
latter point is consistent with experimental observations (ref. 1). 

Poisson's ratio - The frequency and cumulative distribution diagrams for 
the major Poisson's ratio (u^i2) are shown in figure 7 and those for the minor 
(o>q, 21> are shown in figure 8. Examining the diagrams for Case 3 we see that 
the scatter for 1^12 ranges form about 0.35 to 0.75 with a mean value of 
about 0.51. The deterministic value from table II is 0.28 which is less than 
the lower limit in the diagram. The observations to note are: (1) that the 

assumed scatter results in a substantially larger ^12 mean value compared 
to the deterministic value, and (2) that a substantial scatter is probable in 
the measured values for ug.12- The scatter and mean for uq, 21 exhibit simi- 
lar trends as ug,12- This is expected because uq ,21 is not an independent 
quantity; it is calculated from the relationship ^ 2,21 = u &12 1 ^ E C.2 2 - 

Thermal expansion coefficients - The frequency and cumulative distribu- 
tion diagrams for the thermal expansion coefficients are shown in figure 9 for 
longitudinal (<xq,ii), in figure 10 for transverse ( 0 . 2 , 22 ). and ’ n figure 11 for 
the coupling <a&i2>- This last coefficient exists primarily in the presence 
of scatter in the primitive variables which result in an unbalanced equivalent- 
laminate ply. 

It is interesting to note that the diagrams for Case 2 exhibit less scat- 
ter than those for Case 3. This is primarily due to the probable fiber misal- 
ignment. From the Case 3 diagrams, the mean values for a^n, 0222 - and 
a 2,22 are -0.4, 16.5 and 0, pin./in./°F, respectively. The corresponding 
deterministic values are about 0, 18.1 and 0, respectively. The smaller proba- 
bilistic value for ocq, 22 is mainly the result in the scatter of primitive var- 
iables. These diagrams indicate that measured values for <* 0,22 are likely to 
exhibit extensive scatter. 
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It is worth noting that neither the Poisson's ratios nor the thermal 
expansion coefficients for the constituents (primitive variables) are probabi- 
listically defined. Yet the corresponding ply properties are probabilistic 
and with considerable scatter. This demonstrates the complex in situ interac- 
tion and interdependence of the various primitive variables and the need to 
develop a probabilistic approach to composite micromechanics. 


Strengths 

The ply probabilistic strengths are determined by assuming that the ply 
fractured when its weakest subply failed. This is equivalent to first ply 
failure for laminate fracture. 

Longitudinal strengths - The frequency and distribution diagrams for the 
ply longitudinal tensile strength < Sg,i 1 t> are shown in figure 12 and for the 
ply longitudinal compressive strength (Sane* figure 13. As can be seen 
the scatter in Case 3 is: (1) 80 ksi < < 180 ksi with a mean of about 

130 ksi and (2) 40 ksi < Sgjic < 120 ksi with a mean of about 75 ksi. The cor- 
responding deterministic values from table II are 203 and 165 ksi, respec- 
tively. Both of these values are considerably greater than the highest value 
in the scatter range. These results are significant because of the following 
two points: (1) The ply does not fail when its weakest subply fails, other- 

wise the mean would have been close to the deterministic value. (2) The in 
situ mean tensile fracture stress of the fiber is considerably higher than the 
input value. Assuming a proportional (fiber strength/ ply-strength, 203/130) 
the mean for the in situ fiber strength would be about 625 ksi which is about 
56 percent higher than the input mean of 400 ksi. This higher value is repre- 
sentative of very small fiber gage lengths. Both of these lead to the general 
conclusion that ply fracture (strength) is a complex sequence of events which 
generally results in considerable penalty (reduction) relative to in situ fiber 
strength and which is inherently probabilistic. These remarks are applicable 
to the ply longitudinal compressive strength (Sane) as we11 since it exhibits 
similar behavior. 

Transverse strengths - The frequency and cumulative distribution diagrams 
for the transverse tensile ply strength (S&22T) are shown in figure 14 and for 
the compressive i n figure 15. Referring to Case 3 diagrams, it can be 

seen that the scatter is relatively wide: (1) from about 0 ksi to about 7 ksi 

with a mean of 3 ksi for the transverse tensile strength, and (2) from about 0 
to 16 ksi with a mean of 8 ksi for the transverse compressive strength. The 
corresponding values for the deterministic case from table II are 12 and 27 
ksi, respectively. The remarks and discussion made for the longitudinal 
strengths are applicable to the transverse strengths as well, namely: (1) 

transverse ply fracture does not occur when the first subply fails and (2) 
transverse ply fracture is a complex sequence of events which are inherently 
probabilistic. An additional point to be noted is that the deterministic/pro- 
babilistic shift mean ratio is about 4 (12/3) for the transverse tensile ply 
strength and may be indicative of the micro damage tolerance of this 
strength. 

Intralaminar (in-plane) shea r - The frequency and cumulative distribution 
diagrams for the ply intralaminar shear strength (Sg,i2S ) are shown in figure 
16. It can be seen from the Case 3 diagrams that the scatter is relatively 
wide from the 1 to 13 ksi with a mean of about 8 ksi. The corresponding 
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deterministic value from table II is 10 ksi. Again, this ply strength exhibits 
behavior similar to the other four. One noticeable exception is that this ply 
strength has the smallest mean shift ratio 1.25 of all five ply strengths indi- 
cating, perhaps, less brittle behavior for this strength compared to the other 
four. 


Collectively, measured values on any or all of these strengths, would 
exhibit considerable scatter because of their inherent probabilistic nature. 
This scatter would be further accentuated because of the difficulty of perform- 
ing the requisite tests. Though results are not presented here, confidence 
and significance of primitive variables are readily determined using standard 
statistical methods (ref. 3). 


SUMMARY OF RESULTS 

The salient results of an investigation to develop probabilistic compos- 
ite micromechanics are summerized below. 

1. A computational simulation procedure has been developed for probabilis- 
tic composite micromechanics by combining composite mechanics and ply substruc- 
turing, with probabilistic concepts and Monte Carlo simulation. 

2. Anticipated scatter in constituent properties and fabrication processes 
(primitive) variables is directly incorporated in probabilistic composite 

mi cromechani cs . 

3. All the ply mechanical and thermal properties (moduli, thermal expan- 
sion coefficients and strengths) are evaluated simultaneously. This corre- 
sponds to using the same laminate for determining all the properties. 

4. The scatter range and mean in each of the ply properties is presented 
in terms of frequency and cumulative distribution function. 

5. Comparisons of deterministic values with respective probabilistic 
means provide insight into the inherent probabilistic nature of the ply behav- 
ior, anticipated scatter in measured properties, and probable type of fracture. 

6. Ply fracture is a complex sequence of events and is not generally con- 
trolled by the weakest subply and ply fracture stresses will generally exhibit 
wide scatter. 

7. Probabilistic composite micromechanics provide extensive information 
which formally relates ply phenomenological behavior to a large number of com- 
plex and interacting uncertainties at the micromechanics level. 
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TABLE I. - INPUT DATA FOR SAMPLING 


Input 

Determi ni sti c 

Probabi 1 i sti c 


Case 1 

Case 2 

Case 3 

Theta, deg 

0.0 



p 


0.0 

0.0 

o 


5.0 

10.0 

FVR 

P 

0.50 

0.5 

0.5 

o 


0.1 

0.2 

VVR 

\ 

0.01 

0.3 

0.3 

k 


0.03 

0.05 

Ef 1 1 , ksi 

31 000 



p 


31 000 

31 000 

a 


1 500 

3 000 

E f22 . ksi 

P 

2 000 

2 000 

2 000 

a 

ksi 

p 

2 000 

100 



2 000 

2 000 

a 


100 

200 

Gf 2 3> 

P 

1 000 

1 000 

1 000 

O 


50 

100 

'"f 1 2 

M 

0.2 

0.2 

0.2 

a 


0 

0 

“fll. Ppro/°F 
p 

0.2 

0.2 

0.2 

a 


0 

0 

af 22 PP m /°E 
P 

5.6 

5.6 

5.6 

a 


0 

0 

s ft , ksi 
{3 

400 

400 

400 

a 


20 

10 

Sf r, ksi 

400 



(3 


400 

400 

a 


20 

10 

E m . ksi 

P 

500 




500 

500 

a 


25 

50 

P 


0.35 

0.35 

0.35 

a 

— 

0 

0 

“nv Ppm/°F 

36 

36 


P 


36 

a 


0 

0 

S mT , ksi 

15 



13 


1 5 

1 5 

o 


20 

10 

s mc ksi 
13 

35 

35 

35 

a 

— 

20 

10 

S mS , ksi 
3 

13 

13 

13 

a 


20 

10 


10 




Property 


E <?11 

E <?22 

G 612 

2 

y 621 
a 5 1 1 
a £22 
a £l 2 
1 T 

|*nc 

S ^22T 

^ 22 C 

S £12S 












(A) MONTE CARLO SIMULATION OF UNCERTAINTIES IN PRIMITIVE 
VARIABLES. 



(B) QUANTIFICATION OF PLY PROPERTY UNCERAINTIES. 

FIGURE 3. - PROBABILISTIC COMPOSITE MICROMECHANICS 
COMPUTATIONAL PROCEDURE LOGIC DIAGRAMS. 
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(A) CASE 2 HISTOGRAM. (13) CASE 3 HISTOGRAM. 




(C) CASE 2 DISTRIBUTION. (D) CASE 3 DISTRIBUTION. 

FIGURE 6. - SAMPLING RESULTS FOR IN-PLANE SHEAR MODULUS 
(G /12 } * 




(B) CASE 3 HISTOGRAM. 



FIGURE 7. - SAMPLING RESULTS FOR POISSON RATIO (MAJOR 


y £ 1 2 ). 
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.125 .175 .225 .275 .325 .375 0 .2 .4 .6 .8 1.0 

RANGE (E -01) 

(C) CASE 2 DISTRIBUTION. (D) CASE 3 DISTRIBUTION. 


FIGURE 8. - SAMPLING RESULTS FOR POISSON RATIO 
(MIN0R-V£21 ) - 



(A) CASE 2 HISTOGRAM. (B) CASE 3 IIIS10GRAM. 



RANGE (E -0G) 


(C) CASE 2 DISTRIBUTION. (D) CASE 3 DISTRIBUTION. 

FIGURE 9. - SAMPLING RESULTS FOR LONGITUDINAL THERMAL 
EXPANSION (a jR11 ) . 
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RANGE (E 03) 

(C) CASE 2 DISTRIBUTION. (D) CASE 3 DISTRIBUTION. 

FIGURE 12. - SAMPLING RESULTS FOR LONGITUDINAL TENSILE 
STRENGTH <S im >. 





FIGURE 13. - SAMPLING RESULTS FOR LONGITUDINAL COMPRES- 
SIVE STRENGTH CS £11C >. 
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-.02 .02 .06 .10 .14 .18 -.06 -.02 . 02 . 06 .10 .14 


RANGE (E 02) 

(C) CASE 2 DISTRIBUTION. (D) CASE 3 DISTRIBUTION. 

FIGURE 14. - SAMPLING RESULTS FOR TRANSVERSE TENSILE 
STRENGTH (S£ 22T >. 



-.06 .05 .15 .25 .35 .45 -.1 0 .1 .2 .3 .4 


RANGE (E 02) 

(C) CASE 2 DISTRIBUTION. (D) CASE 3 DISTRIBUTION. 

FIGURE 15. - SAMPLING RESULTS FOR TRANSVERSE COMPRESSIVE 
STRENGTH (S^O. 
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